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ABSTRACT 



£> ' In this paper we report on the broadband X-ray properties of a complete sample of 

\o : 

t^J- , absorbed Seyfert galaxies hard X-ray selected with INTEGRAL. Our sample is com- 



posed of 33 sources, of which 15 are newly discovered AGN above 20 keV (IGR 
sources) while 18 are already known type 2 AGN ("known")- For 17 sources (15 IGR 
+ 2 "known" sources) we have performed a broadband analysis using both XMM- 
Newton, and INTEGRAL-IBIS data. To have a full view of the complete sample we 
have then complemented the analysis of the 16 remaining sources with already ex- 
isting broadband studies in the same range. The high quality broadband spectra are 
well reproduced with an absorbed primary emission with a high energy cutoff and 
its scattered fraction below 2-3 keV, plus the Compton reflection features (Compton 
hump and Fe line emission). This study permitted a very good characterization of the 
primary continuum and, in turn, of all the spectral features. 

A high energy cut-off is found in 30% of the sample, with an average value 
below 150 keV, suggesting that this feature has to be present in the X-ray spectra 
of obscured AGN. The hard X-ray selection favours the detection of more obscured 
sources, with the log N# average value of 23.15 (standard deviation of 0.89). The 
diagnostic plot N# vs F eorr (2-10 keV)/F(20-100 keV) allowed the isolation of the 
Compton thick objects, and may represent a useful tool for future hard X-ray obser- 
vations of newly discovered AGN. We are unable to associate the reflection compo- 
nents (both continuum and Fe line) with the absorbing gas as a torus (as envisaged in 
the Unified Model), a more complex scenario being necessary. In the Compton thin 
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sources, a fraction (but not all) of the Fe K line needs to be produced in a gas located 
closer to the black hole than the Compton thick torus, and this is possibly associ- 
ated with the optical Broad Line Region, responsible also for the absorption. We still 
need a Compton thick medium (not intercepting the line of sight) likely associated 
to a torus, which contributes to the Fe line intensity and produces the observed re- 
flection continuum above 10 keV. The so-called Iwasawa-Taniguchi effect can not be 
confirmed with our data. Finally, the comparison with a sample of unobscured AGN 
shows that, type 1 and type 2 (once corrected for absorption) Seyfert are charac- 
terized by the same nuclear/accretion properties (luminosity, bolometric luminosity, 
Eddington ratio), supporting the "unified" view. 

Key words: galaxies: active - galaxies: Seyfert - X-rays: galaxies 

1 INTRODUCTION 

According to the widely accepted Unification Model (UM) for active galactic nuclei ( Antonucci 
1993, Urry & Padovani 1995), type 2 and type 1 AGN are believed to share the same central region 
and emission mechanisms with the only difference being due to orientation. In Seyfert 2, only 
optical narrow emission lines are visible, while in Seyfert 1 both narrow and broad emission lines 
are detected (Antonucci & Miller 1995). Type 2 AGN are seen edge-on, across a large amount of 
obscuring gas that prevents us from directly observing the central core and the broad line regions 
(BLR). This scenario is strongly supported by spectropolarimetric observations of hidden broad- 
line regions (HBLRs) in several type 2 AGN, and by the X-ray observations, demonstrating that 
Seyfert 2 galaxies usually have absorption columns greatly exceeding the Galactic ones (Cappi et 
al. 2006). Despite this, and in view of several pieces of evidence coming from recent observations 
of different classes of AGN, the UM is still being debated. In particular, the association of the 
Compton thick (a cold gas with N H > cr~j - 10 24 cm 2 ) obscuring medium with a uniform torus 
at intermediate scale between the BLR and the narrow line regions (NLR), has been questioned by 
the presence of column density variability on timescales of hours, requiring an X-ray absorber no 
larger than the BLR. Such changes have been detected in the brightest absorbed Seyfert galaxies, 
NGC 1365, NGC 4151 and NGC 4388 (Risaliti et al. 2007, Puccetti et al. 2007, Elvis et al. 2004). It 
has also been proposed that the co-existence of a Compton thick torus and Compton thin material, 
extended on a much larger scale, better accounts for the observed phenomenology in absorbed 
AGN (e.g. Matt 2000). In addition, some Compton thick objects exhibit no optical activity and a 

* E-mail: alessandra.derosa@iasf-roma.inaf.it 
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weak scattering X-ray component, indicating that the torus is geometrically thicker than previously 
thought (Ueda et al. 2007). These absorbing media are also responsible for the reflected features 
(lines and continuum) observed in the X-ray energy domain: a Compton hump above 10 keV 
and a Fe line at 6.4-6.9 keV (Perola et al. 2002). Also, the physical/geometrical properties of 
environmental materials are subjected to the luminosity of the central supermassive black-hole 
(SMBH) (Bianchi et al. 2007). On the other hand, it is also important to understand the behaviour 
of the complex multi-component spectrum of an AGN to clarify the physical view of the central 
engine. To this end, correct modelling of the X-ray continuum emission is mandatory. Finally, the 
study of obscured AGN represents a powerful tool to trace the growth of SMBH and to understand 
how accretion took place throughout cosmic history. In particular, Compton thick AGN are the best 
candidates to explain the unresolved fraction of the Cosmic X-ray Background (CXB) in the X-ray 
energy domain around 30 keV (Gilli et al. 2007). In this respect, hard X-ray observations, above 
10 keV, represent a major improvement in selecting obscured AGN that are biased by absorption 
much less severely than those selected in the 2-10 keV band. However, because of the limited 
sensitivity of current X-ray detectors at such high energies, these studies can only be conducted in 
the very local Universe. 

This kind of study has been performed in a statistical way with a variety of different samples 
over a broadband energy range: BeppoSAX (Risaliti 2002), INTEGRAL and XMM-Newton (De 
Rosa et al., 2008), Suzaku and Swift-BAT (both type 1 and type 2 AGN, Eguchi et al. 2009, 
Winter et al. 2009, Fukazawa et al. 2011). In the hard X-ray energy range, the IBIS instrument 
(Ubertini et al. 2003) on board INTEGRAL (Winkler et al. 2003) and the Burst Alert Telescope 
(BAT, Barthelmy et al. 2005) on board Swift (Gehrels et al. 2004) are currently performing surveys 
of the hard X-ray sky, discovering large numbers of sources seen for the first time at these energies 
(Malizia et al. 2009, Landi et al. 2010, Tueller et al. 2010, Winter et al. 2009, Bassani et al. 2006). 
The 4th INTEGRAL/IBIS catalogue (Bird et al. 2010) includes more than 700 hard X-ray sources 
detected over a large part of the skjQ above 20 keV and with positional accuracies ranging between 
0.2 and 5.5 arcmin, depending on the source significance. 

To obtain the best quality broadband high energy data, the INTEGRAL/ AGN team has under- 
taken a programme of follow-up observations of hard X-ray selected AGN with XMM-Newton 
(Molina et al. 2008, Panessa et al. 201 1) 

1 70% of the INTEGRAL hard X-ray sky around the Galactic plane region is covered to better than 1 mCrab sensitivity (~ 2xl0~"erg cm~ 2 s ), 
compared to the 5 mCrab level of 90% of the extragalactic regions (see details in Bird et al. 2010). 
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In this paper we report on the broadband spectral properties of all type 2 AGN included in 
the hard X-ray selected AGN complete sample as defined by Malizia et al. (2009). This paper is 
the continuation of the work we presented for a smaller sample of 7 hard X-ray selected type 2 
AGN (De Rosa et al. 2008). The actual sample comprises 33 AGN optically classified as type 
1.8-2. The new 15 sources detected for the first time above 20 keV with INTEGRAL have been 
observed within our XMM-Newton follow-up campaign, while we make use of published data to 
characterise the broadband spectra of 18 previously known objects. 

The main goals of this study are (1) to measure the spectral properties of obscured AGN: the 
intrinsic continuum shape (both photon index and high energy cut-off) together with the absorption 
distribution and (2) to accurately constrain the strength of the reflection components: the Fe K line 
and the Compton hump. The line properties tell us about the obscuring/reflecting medium while 
the Compton continuum is an important yet uncertain parameter (like the high energy cut-off) in 
population synthesis models of the CXB (Ueda et al. 2003, Gilli et al. 2007). 

The sample is defined in Sect. [21 details of the INTEGRAL and XMM-Newton observations 
are given in Sect. [31 while the broadband spectral analysis is described in Sect. 0] The results are 
presented and discussed in Sect. [51 while in Sect. [6] we investigate the accretion properties of the 
sample compared to a sample of type 1 hard X-ray selected AGN. Our results are summarized in 
SectE 



2 THE SAMPLE 

From the list of INTEGRAL selected AGN included in the 3rd IBIS survey (Bird et al. 2007), a 
complete sample has been extracted by means of the V/V max test, which was first introduced by 
Schmidt (1968) as a test of uniformity of the distribution in space for a flux-limited sample of 
objects (see Malizia et al. 2009 for details). There are 88 objects detected in the 20-40 keV band 
forming this complete sample of INTEGRAL selected AGN: 46 objects are of type 1 (Seyfert 1- 
1.5, of which 5 are Narrow Line Seyfert Is) and 33 of type 2 (Seyfert 1.8-2); only 9 Blazars (BL 
Lac-QSO) are included in the catalogue. The list of objects we discuss in this paper is composed 
of all the type 1.8-2 objects in the INTEGRAL complete sample of AGN defined in Malizia et 
al. (2009) except for NGC 1275 and NGC 5506: the first source has been excluded because of 
the Perseus cluster contamination in the low energy band, while the second source has been re- 
classified as a Narrow Line Seyfert 1 (see Panessa et al. 2011). To the present sample we have 
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added two sources (IGR J00040+7020 and IGR J20286+2544) because they have been observed 
with moderately long exposure with XMM-Newton and are hard X-ray INTEGRAL selected. 

All sources are optically classified, and their redshift distribution has an average value of 0.03 
( z m«^=0.252). Within our complete sample, 15 are IGR sources, i.e. AGN first detected by IN- 
TEGRAL above 20 keV, while 18 are already known type 2 Seyfert (hereafter "known"). For the 
newly discovered objects we started an XMM-Newton Large Programme (PI. L. Bassani) with 
the primary goal to characterize their broadband spectra. The last XMM-Newton observation was 
performed in February 2011. For the known AGN we used results of broadband studies already 
published, when available, with Suzaku, INTEGRAL and XMM-Newton, INTEGRAL and Chandra 
or BeppoSAX. When broadband data was completely missing from the literature, we analysed 
unpublished XMM-Newton and INTEGRAL data. In conclusion, for this work we performed data 
analysis and a spectral study of all IGR sources plus 2 known objects that had no broadband study 
already published (NGC4507 and NGC7172) with the XMM-Newton data. The sample, together 
with additional information (observation dates, exposures in the low and high energy domains, 
2-10 and 20-100 keV counts/s and, for the known sources, the instrument used with relative ref- 
erence), is presented in Table \T\ where the "new" IGR sources are described in the upper part, and 
the already "known" sources in the bottom part. 

3 OBSERVATIONS AND DATA REDUCTION 
3.1 INTEGRAL/IBIS 

Though the sources reported here have been hard X-ray detected in the 3rd INTEGRAL survey 
(Bird et al. 2007), to better characterize the spectrum above 20 keV we use the data available in 
the 4th INTEGRAL catalogue (Bird et al. 2010), making use of the longer on-source exposures. 
We also stress here that we used INTEGRAL/IBIS data processed by the survey team following the 
procedures as described below and reported in Bird et al. (2010). For each source, they combined 
INTEGRAL/IBIS data from several pointings performed between revolution 12 and 530 to pro- 
vide sufficient statistics for weak mCrab flux sources. First they generated IBIS/ISGRI images for 
each available pointing in 13 energy bands with the ISDC offline scientific analysis software OS A 
version 7.1 and they then extracted count rates at each source position from individual images to 
provide light curves in the various energy bands sampled; since the light curves did not show any 
sign of variability or flaring activity, average fluxes were then extracted in each band and combined 
to produce the source spectrum. In Table[H we report the details of the IBIS/ISGRI observations; 
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Table 1. Journal of the observations 



Name 


z 


*(RA,DEC) 


Date 


Exp Time 


Cts(2-l()keV) 


Cts(20-100 keV) 








XMM obs 


XMM-pn/IBIS 






IGR J 




(J2000) 




(ks) 






00040+7020 


0.096 


00 04 01.93 ;70 19 18.6 


2008-12-29 


18/2614 


0.326±0.004 


0.23±0.02 


m^A^ ZL 1 1Q 
\J 1 JOJ-41 Jy 


n no 1 


\) 1 JO iy.04 ,-41 3 1 4Z.1 


zUUo-l l-lo 


17/1 ^SO 


n 'ios.-i.n nn^ 

U. Jzo±U.UUj 


U.ZU±U.Uj 


09523-6231 


0.252 


09 52 17.00 ;-62 31 00.0 


2009-01-20 


17/1607 


0.321 ±0.005 


0.20±0.03 


1009-4250 


0.033 


09 52 17.00 ;-62 31 00.0 


2007-06-14 


21/1020 


1.936±0.003 


0.40±0.04 


1 f \ A A A 

10404-4625 


/"\ a o o a 

0.0239 


1 A A A O O O 1 A £1 O C O /~ C 

10 40 22.31;-46 25 26.5 


o AA/ ii OA 

2006-1 1-29 


1 1/781 


1 O O i A A 1 

1.23±0.01 


A A O i A A C 

0.48±0.05 


1 iai r co An 

12026-5349 


0.03 


1 O AO in /O CO CA AC n 

12 02 47.68 ;-53 50 08.3 


OAAA 1 O OA 

2009-12-30 


O C I A O A C 

25/1805 


A ^1 /~ , A AAC 

0.726±0.005 


A /" A i A AO 

0.60±0.03 


1 o aa 1 . i in/j nni \ 

13091 + 1137(4992) 


0.03 


1 O A A AC C O 1 1 O O AO 1 

13 09 05.58 ;11 38 02.7 


OA A/""" A/""" 0*"7 

2006-06-27 


13/495 


A 1 n/** i A AAO 

0.106±0.003 


A C 1 i A A/~" 

0.51 ±0.06 


\ A C A C C C ,1 O 

14515-5542 


0.018 


1 A ci o o a f cc An no a 

14 51 33.16; -55 40 38.8 


■OA 1 A AO A 1 

2010-02-01 


18/2617 


A CCC . A AA/T 

0.555±0.006 


A OO . A AO 

0.28±0.03 


16024-6107 


0.011 


1 /" Al -lOOO C 1 A O C ,1 

16 01 48.23; -61 08 54.7 


OAAA A A A A 

2009-09-09 


1 A 'OAOO 

19/2022 


A O CO . A AAC 

0.358±0.005 


A OA . A AO 

0.20±0.03 


1 635 1 -5 806( 1 3 /-(j34) 


0.009 1 


1/ 1ft iA . rti (\r r\A 1 1 

16 35 10.40 ;-58 06 04.0 


2006-02- 1 3 


1 //2351 


n mn . A aa^ 

0.039±0.002 


A T7 i A AQ 
0.2/±0.03 


1 ""7C 1 O 1A1 1 

17513-2011 


0.047 


17 51 13.62 ;-20 12 14.6 


TAA/" AA A A 

2006-09-04 


O /CAOA 

8/5980 


A C <"7 . A AAA 

0.567±0.009 


A A O . A AO 

0.42±0.02 


i An An 

20186+4043 


0.0144 


OA 1 o OO CA An A 1 AA I"1 

20 18 38.50 ;40 41 00.0 


OA1A 11 A1 

2010-1 1-01 


1 O IA OAO 

12/1202 


A An A .A AA/T 

0.434±0.006 


A OA . A AO 

0.29±0.03 


TATn X . n C A A 

20286+2544 


0.01 


OA OO OC AC OC /I /I AA "1 

20 28 35.06 ;25 44 00.2 


OAA/""" A A OO 

2006-04-23 


1 1/469 


A 1 A O . A AA A 

0.142±0.004 


A CO i A A A 

0.53±0.04 


oooao i nion 

23308+7120 


0.037 


O O O A O ^7 O C 11 O O /I ,1 O 

23 30 37.25 ;71 22 44.8 


OA 11 A1 OO 

2011-01-23 


O O / 1 A A C 

33/1995 


A 1 S O . A AAO 

0.162±0.002 


A A A i A AO 

0. 14±0.03 


O O CO ,1 i C O ,1 O 

23524+5842 


0.164 


T) CI 0"7 AA CO 1 O AA A 

23 52 27.00 ;58 42 00.0 


OA1 1 AO OA 

20 l 1-02-20 


1 o /o a net 

18/3498 


A OOO i A A A A 

0.232±0.004 


A 1 ^7 i A AO 

0. 17±0.02 


known 






*Date 




^ Instrument 


Reference 


NGC788 


0.0136 


02 01 06.45,;-06 48 57.0 


1999-01-16 


39/594 


Chandra - INTEGRAL 


De Rosa et al. 2008 


NGC1068 


0.038 


02 42 40.77 ; -00 00 47.8 


1996-12-30 


95/65 


Beppo-SAX 


Matt et al. 1997 


NGC1142 


0.0288 


02 55 12.27; -00 11 01.8 


2007-01-23 


102/81 


Suzaku 


Eguchi et al. 2009 


Mkn3 


0.0135 


06 15 36.46 ; 71 02 15.2 


2005-10-23 


87/81 


Suzaku 


Awaki et al. 2008 


MCG-05-23-16 


0.0085 


09 47 40.17; -30 56 55.9 


2001-12-07 


98/71 


Suzaku 


Reeves et al. 2007 


NGC3281 


0.0115 


10 31 52.09 ; -34 51 13.4 


2000-05-20 


71/35 


Beppo-SAX 


Vignali & Comastri 2002 


NGC4388 


0.0084 


12 25 46.82; 12 39 43.4 


2001-06/07 


13/512 


XMM-Newton- INTEGRAL 


Beckmann et al. 2004 


NGC4507 


0.0118 


12 35 36.62 ; -39 54 33.2 


2001-04-04 


36/393 


XMM-Newton- INTEGRAL 


this work 


LEDA170194 


0.036 


12 39 06.30; -16 10 47.2 


2005-07-25 


3/200 


Chandra -INTEGRAL 


De Rosa et al. 2008 


NGC4945 


0.0019 


13 05 27.28 ; -49 28 04.4 


2006-01-15 


99/69 


Suzaku 


Itoh et al. 2008 


CenA 


0.0018 


13 25 27.62 ; -43 01 08.8 


2005-08-19 


70/60 


Suzaku 


Markowitz et al. 2007 


NGC5252 


0.023 


13 38 15.87; 04 32 33.0 


1998-01-20 


61/29 


Beppo-SAX 


Risaliti 2002 


Circinus 


0.0014 


14 13 09.91 ; -65 20 20.5 


1998-03-24 


137/63 


Beppo-SAX 


Mattetal. 1999 


IC4518A 


0.0163 


14 57 40.50 ; -43 07 54.0 


2006-08-07 


11.4/898 


Chandra - INTEGRAL 


De Rosa et al. 2008 


NGC6300 


0.0037 


17 16 59.47 ; -62 49 14.0 


2001-1999 


44/77 


XMM-Newton- Beppo-SAX/PDS 


Matsumoto et al. 2004 


ESO103-G35 


0.0133 


18 38 20.36; -65 25 39.2 


2002-03-15 


12/44 


Chandra- INTEGRAL 


De Rosa et al. 2008 


CygA 


0.0561 


19 59 28.36; 40 44 02.1 


2000-05-21 


9/29 


Chandra - RXTE/PCA 


Young et al. 2002 


NGC7172 


0.0087 


22 02 01.90 ; -31 52 11.6 


2007-04-24 


28/250 


XMM-Newton- INTEGRAL 


this work 



* Optical position in hh mm ss.ss; dd mm ss.s. tFor Suzaku observation exposures refer to XIXS and PIN data; for Beppo-SAX observation 
exposures refer to MECS and PDS data; for Chandra observations exposures refer to ACIS data; for XMM-Newton observation refer to pn data; 
for RXTE to PCA data. ? Date of observations performed with INTEGRAL and XMM-Newton or Chandra refer to XMM-Newton/Chandra 

observations. 

(optical) best fit positions, exposures and count rates in the range 20-100 keV band. 



3.2 XMM-Newton 

To avoid cross-calibration problems and for the sake of simplicity, we discuss in this paper only 
the data of the EPIC/pn camera (Struder et al. 2001), however the MOS1 and MOS2 spectra 
have been checked for consistency. Data were reduced using SAS vlO.O employing eproc with 
standard settings and the most updated calibration files available at the time of the data reduction 
(April 201 1). Source spectra were extracted from circular regions surrounding the source position 
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while background spectra were extracted in the same CCD chip from circular regions free from 
contaminating serendipitous sources. Each spectrum was rebinned in order to have at least 30 cts 
per bin allowing us to apply x 2 statistics. Spectral fits were performed in the 0.3-10 keV energy 
band. The information on the XMM-Newton observations (date, exposure and 2-10 keV count 
rate) is shown in Tabled] 



4 SPECTRAL ANALYSIS 

We fitted simultaneously the soft and hard X-ray spectra available with INTEGRAL/IBIS and 
XMM-Newton in the 0.3-100 keV energy range. To take into account possible flux variation of 
the source between the non simultaneous XMM-Newton and INTEGRAL data, we left a cross- 
calibration constant, between IBIS and pn, free to vary during the fit procedure. Possible miscal- 
ibration between INTEGRAL and XMM-Newton could mimic or hide the presence of a Compton 
reflection component above 10 keV; we will discuss this point in Sect. 15 .31 

For the sake of clarity and to ease the comparison with the published results on the known 
objects, we considered a baseline model (BLM thereafter) which is commonly used to fit the 
broadband X-ray spectra (0.3-100 ekV) of absorbed AGN (Comastri et al. 2010, Ueda et al. 2007, 
Risaliti 2002, De Rosa et al. 2008). The BLM is composed of: 1) a primary continuum extending 
up to the hard X-ray energy domain modelled by a cutoff power-law having photon index Y and 
high energy roll-over E c ; 2) a soft X-ray scattered component modelled by a power-law having 
photon index Y s assumed to be equal to Y and a scattering fraction f sc ; 3) a Compton reflection 
continuum above 10 keV with the inclination angle 8 of the reflector fixed to 45 deg, and solar 
abundance, plus the Fe Ka emission line that is reproduced with a gaussian profile. We assume 
the Fe line to be intrinsically narrow (<x <0.01 keV, FWHM<1000 km s _1 ) and no additional 
Fe K/? emission component is included in this model; this is due mainly to the fact that we are 
characterizing the continuum properties. However, in some cases we need to relax the assumption 
about narrowness to well reproduce the Fe line; these cases will be discussed in Sect. 15.41 All the 
spectral components are absorbed by the Galactic column density along the line of sight N^° 
while the primary power-law is also absorbed by a second cold medium with column density N H 
We can write the BLM of the photon spectrum F(E) as follows: 

e -<r th N%">( e -cr tH N HKE -r e -E/E c + f sc KE~ r > + R(E) + G(E) + S(E)) 



1 



2 ]\jGn/ va [ ues have calculated through the web interface http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl using the Dickey & Lockman 
(1990) HI in the Galaxy. 
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where G(E) represents the Fe-K Gaussian emission line, R(E) is the Compton reflection com- 
ponent (with the shape fixed to that of the primary emission), and S(E) represents any additional 
soft components needed to reproduce the soft X-ray spectra (see Sect. 15.51) . The total baseline 
model (excluding Galactic absorption) is expressed as zphabs*cutoffpl + const*cutoffpl 
+ zgauss + pexrav + soft comp in the XSPEC terminology. The broadband spectra of IGR 
sources are shown in Fig|9l \W\ and QT| while the best fit parameters for the whole sample are 
reported in Tabled 

In this paper we quote statistical errors at the 90% confidence level for one interesting parame- 
ter unless otherwise specified. In order to calculate the luminosities we use a standard A cold dark 
matter cosmology with the following cosmological parameters: H =70 km s" 1 Mpc~ ! , A =0.73, 
a M =0.21 (Bennett et al. 2003). 

5 RESULTS AND DISCUSSION 
5.1 The primary continuum 

The BLM reproduces quite well all the spectra analysed since the ^f 2 /dof ranges between 0.77- 
1.12. In Fig. \T\ we show the distribution of the photon index of the primary absorbed power-law 
for the whole sample (on the left) as well as for the IGR and known objects separately (on the right) 
to search for possible differences in the radiative emission mechanisms between the two sample 
of sources. The mean value of the photon index and its standard deviation for the whole sample 
is r=1.68, cr=0.30. This value is in very good agreement with those found in previous analyses 
performed in a limited energy range (2-10 keV) in both type 1 and type 2 AGN (e.g, Singh et 
al. 2011, Bianchi et al. 2009). A detailed comparison between type 1 and type 2 Seyfert will be 
presented in Sect.[6l We note that the average value of the IGR population is only marginally flatter 
than the known one, being r /GR =1.63, cr /GR =0.34 and T known = l.lA, cr known =0.2l. A Kolmogorov- 
Smirnov test applied on these two populations indicates that they come from the same parent 
distribution with a probability of P=0.385. To explore a possible correlation between the primary 
continuum and the brightness, we plot in Fig. [2](left panel) the photon index against the 2-10 keV 
luminosity corrected for absorption: we do not find any evidence of correlation, a Pearson test for 
linear fit yields r=-0.02. 

We measured the high energy cutoff in 10 out of 33 sources (i.e. 30 per cent of the sample), 
namely: NGC4507, IGR J 12026-5349, IGR J13091 + 1 137, IGR J201 86+4043, IGR J23308+7120, 
IC4518A, NGC7172, Circinus, NGC788 and ESO103-G35, while only an upper limit at 300 keV 
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Table 2. Broadband baseline model (BLM) composed by an absorbed power-law with exponential cutoff+reflection continuum+ Fe line + soft 
scattered component with r so f t =T 



Name 
IGRJ 



1 




3 R 




5 f 
•sc 


6cw 
IJ. w 


1 44+O.O6 

-0 02 


33+J 


< 1.8 


> 120 


9+° 5 

-0 5 


50+ 48 

-47 


| 7 g+o:08 
* -0.03 


7 9+0.4 
'■ z -0.9 


■X c+0.6 


> 100 


5+ 4 ' 

—2 


163+' 13 

1 -113 


i 7 c +0.09 
-0.09 


63+. 4 


<2 


>40 


1 8+°-3 


77+35 


1 66+ - 04 


257 +15 


1 4+ 03 

1- -0.3 


> 150 


1 0+o 1 

-0.1 


70+| 


9 17+0.03 
-0.03 


41 + [ 


4 3+ L4 


> 250 


<0.2 




1 16+ 005 

1.1D _ 003 


29% 


1 1+ - 5 
• -0.5 


76+ 14 
'°-ii 


4 l+o 1 


100 -59 


i i 5 +0.09 


439+36 
-,3fi 


< 0.4 


10 7+193 


0.5+ 01 


379+81 


164 +8:°? 


3 3+'* 


1 4+ L9 


> 80 


30 + ^° 


182+| 


9 26 +0.08 
zzo -0.10 


2 5+ 03 


J - J -0.6 


> 70 


< 3.2 


121+ 80 


j 2+0.1 


2000 + |0 


0.8+0-2 


> 150 


9 0+o 6 


1400+ 300 


1 

6+8:1 

'"-0.2 


58+5 


3.5+ 3:8 

2 0+' :3 ' 
^""-0.9 


> 150 


< 12 
1 Q+0.9 
' -0.4 


< 200 
80+so 


1 4+0.1 


574^ 


< 1 


> 180 


5+°- 2 


112+1 


2 +0 - 2 




< 2.7 


63 +340 


< 2 


100+ 


1 5+ 0J 


2 9 +4 


<3 


>120 


5+2.0 
-2.0 


100+* 

— OU 



7 F2-10keV 8 F20-100keV 



' log L ( 20-100faV) 



00040+7020 
07565-4139 
09523-6231 
1009-4250 
10404-4625 
12026-5349 
13091+1137 
14515-5542 
16024-6107 
16351-5806 
17513-2011 
20186+4043 
20286+2544 
23308+7120 
23524+5842 



0.35 
0.32 
0.37 
0.31 
1.2 
0.85 
0.21 
0.53 
0.30 
0.05 
0.58 
0.53 
0.28 
0.14 
0.28 



1.6 
1.3 
1.3 
2.7 
2.7 
3.7 
3.8 
1.9 
1.3 
1.5 
3.0 
2.0 
5.5 
0.9 
1.1 



44 
42.5 

45 
43.3 
43.3 
43.3 
43.1 
42.6 
41.9 
42.1 
43.5 
42.6 
42.8 
42.9 
44.4 



44.6 
43.1 
45.5 
43.9 
43.6 
43.8 
43.8 
43.2 
42.6 
42.5 
44.2 
42.9 
43.3 
43.4 
44.9 



known 



NGC788 
NGC1068 
NGC1142 
Mkn3 
MCG-05-23-16 
NGC3281 
NGC4388 
NGC4507 
LEDA170194 
NGC4945 

CenA 
NGC5252 
Circinus 
IC4518A 
NGC6300 
ESO103-G35 
CygA 
NGC7172 



■ 9<;+0.05 
i ""-0.17 
2 13+0- 17 
1 7g+8:o? 

1 00+8:86 
' -fl-RS 

1 95+0.03 
-0.03 
1 95+O.I8 
1? -0.18 

i.8or™ 
1.7* 

1.6+°' 
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1 oq+8:°3 
l.S3_ p 
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, 5+-8J 1 

l.J _ 03 
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j 9+0.2 

1 c 0+8:0? 

1.30 An , 
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le4 
631+. 2 ' 
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16.5+ - 3 



1500 
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-0.3 
+200 
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53 00 +400 
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4300+ 4 oo 

140! 
240+'°, 

189+* 

200+ 10 

-jo 

8+0.7 
-0.7 



9+ 13 
u.5»_ 7 



j 49+0.15 
1 36+o^ D 

• -0.2 

3+°- 1 

—0.1 

3+°- 2 
1.5 
<4.8 



< 0.2 

1 6+ L2 

-1.2 

2+0- 1 

2 6+ 3 - 5 

^.u_ 2 - 

1.1+E 



1.7 



73.: 



-0.5 



3+0- 1 

u "'-0.1 



62+ 38 

DZ -24 



> 300 

> 200 

> 170 

> 60 

> 180 
126+ 153 

>210 

> 80 
>400 

> 50 
56+ 8 

-23 
7Q+6O 
' -30 

>250 
50+ 250 



66+ 12 



<2 

26+0 05 
-955 

0.9+0; 

50+^ 01 
u.ju _ 001 



1.5 



<1.6 



0.2 



1.1+0.1 

20+^0' 

u - zu -0.03 



5 39+ 231 
J -109 

1000+ 140 


0.61 
0.52 


3.8 
3.9 


42.8 
44 


43.1 
41.9 


218+ 


0.42 


5.3 


43.7 


44.0 


855+ 32 


0.64 


10 


43.2 


43.6 


70+° 


8.76 


14.5 


43.2 


43.5 


1180! 4 oo 
29Q+10 

UT® 

11 -27 
< 3000 


0.29 
2.3 
1.3 
2.0 


7.2 
20 
17.4 
4.9 


43.2 
43 
43.2 
43.9 


43.3 
43.6 
43.8 
44.2 


1600+ 300 
lu -300 

83+| 


0.54 
21.2 


28.9 
64 


43 
41.9 


42.6 
43.0 


278+ 2 '5 

2250+10 
ZZJU -100 

360+ 106 
140+50 
1 -50 
395+100 
J -100 

182+ 43 


0.26 
1.4 
0.29 
0.86 
2.4 
1.2 


0.66 
23 
2.9 
9.9 
7.9 
4.0 


42.7 
42.2 
42.6 
41.8 
43.4 
44.6 


42.9 
42.3 
43.3 
42.1 
43.5 
44.5 


76+ '0 


4.3 


8.2 


43 


43.0 



* frozen in the fit; 1 photon index of the primary absorbed, power-law; ; 2 absorber column density in 10 21 cm" 2 ; 3 Relative reflection; 4 high 
energy cutoff in keV; 5 scattered fraction 
in the 2-10 keV energy range in 10~" 



6 Fe Ka equivalent width measured in eV with respect to the observed total continuum; 7 Observed flux 

10~ 

2-10 keV energy range corrected for absorption in erg s ; 



erg cm 2 s 1 ; 8 Observed flux in 20-100 keV energy range in 10 1 1 erg cm 2 s 1 ; 9 Log of the luminosity in 



has been found in IGR J 1635 1-5806. In the remaining sources we obtained lower limits to E c . 
The average value of the ten measured high energy cutoffs is E c =73 keV with a standard deviation 
cr=25 keV. In Fig. [2] (right panel) we plot the high energy cutoffs as a function of the photon index 
T. All the measured values are below 150 keV, while lower limits (filled triangles) are well grouped 
below 300 keV. These values are consistent with those observed in the Seyfert 1 by BeppoSAX 
(Perola et al. 2002, De Rosa et al. 2007). We do not find any evidence of correlation between Y 
and E c , the correlation coefficient being r=0.1. Our analysis strongly suggests that the high energy 
rollover is an ubiquitous property of Seyfert galaxies, and that its value ranges in the same intervals 
for type 2 and type 1 Seyferts (Molina et al. 2008). It is worth noting that we were able to measure 
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Figure 1. Left panel: Primary photon index T distributions for the whole sample. Right panel: Photon index distributions for two separated popula- 
tion: IGR (diagonal dashes) and known literature sources (horizontal dashes). 
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Figure 2. Left panel: photon index for the whole sample as a function of the 2-10 keV luminosity corrected for absorption. Filled boxes are IGR 
sources and open circles are known sources. Right panel: High-energy cutoff vs photon index for the whole sample. Filled triangles are lower limits 
to E c while the filled dot is the upper limit of IGR J16351-5806. 



E c in the objects with the highest S/N ratio, suggesting that future deep hard X-ray observations 
will be the right way to investigate the presence of this component in the spectra of AGN. 



5.2 The absorption properties 

The absorbing column density distribution of the primary emission is shown in the left panel of 
Figf3l Its average value is 23.15 with a standard deviation <x=0.89. Excluding the six Compton 
thick sources (IGR J16351-5806, NGC 3281, NGC 4945, Circinus, NGC 6300, NGC 1068) we 
obtain an average value (standard deviation) of 22.85 (0.670). The absorbing column density dis- 
tribution of the whole sample is in full agreement with that reported by Malizia et al. (2009) for 
the same sample. When our N H distribution (left panel in Fig. [3]) is compared with the absorption 
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Figure 3. Left panel: Distribution of the absorbing column density of the whole sample. Right panel: The column density distribution as a function 
of the softness ratio F OJS (2-10 keV)/F(20-100 keV). The lines represent a region in the diagnostic defined when accounting for the typical values 
of the photon index and high energy cutoff, with errors, as found in our analysis. 

A clear trend is evident with Compton thick sources showing a lower softness ratio. 



distribution obtained with 2-10 keV data only, e.g. in the XMM-Newton and Chandra deep surveys 
(Giacconi et al. 2002, Worsley et al. 2004, Akylas et al. 2006), it becomes evident that the hard 
X-ray data strongly favours the detection of more heavily obscured sources than 2-10 keV data 
only. Nonetheless, in the hard X-ray domain above 10 keV a significant fraction of the source flux 
might be lost if the source is heavily Compton thick (log(N//) >25). In fact, a number of objects 
of this type is not present in our sample as well as in the Swift/BAT survey above 10 keV (Burlon 
etal. 2011). 

In Fig. [3l right panel, we plot the absorbing column density as a function of the ratio of the 
observed 2-10 keV flux versus the 20-100 keV flux. This diagnostic plot was proposed by Malizia 
et al. (2007) in an attempt to isolate peculiar objects among a sample of 34 INTEGRAL AGN 
observed with Swift-XRT. The plot shows a clear anticorrelation between N H and the softness 
ratio: this ratio becomes lower as the column density increases. The lines define a region in the 
plot obtained by assuming an absorbed power-law with high energy cutoff, where we have used the 
average values of T and E c , with errors, as found in our sample (see Sect. 15.11) . In this diagnostic 
plot, Compton thick AGN are well isolated, all six Compton thick objects of our sample are located 
above the line. 



5.3 The reflection continuum 

Our BLM includes the Compton reflection components, namely, the Compton hump and the Fe 
Ka line. We have measured the reflection fraction (R) in 22 sources (average value R=1.63 with 
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cr=1.15), while for 8 objects we only found upper limits. The distribution of the values of R derived 
for the whole sample is shown in Fig. 0] where upper limits are shown as obliques bars while the 
measured values are filled with horizontal bars. 

In the right panel of Fig. |4] we show the relative reflection R as a function of the photon index 
for the whole sample. The data show a poor correlation coefficient (r=0.4), indicating that the 
correlation claimed by Zdziarski et al (1999) between R and Y is not found in our data. It is worth 
noting that the fit of our BLM produces a degeneracy between the reflection fraction R and the 
photon index Y: a higher Y produces a higher value of R. Another point to stress is that the soft 
and high energy data (below and above 10 keV) are non-simultaneous. INTEGRAL/IBIS data are 
the result of the average flux state of the source over a time interval of about 5 years (the time 
interval covered by the fourth INTEGRAL catalogue Bird et al. 2010). XMM-Newton observations 
last at most 20 ks and are therefore a snapshot not necessarily simultaneous with INTEGRAL/YBIS 
pointings (see Tabled)). To take into account possible flux variation between XMM-Newton and 
INTEGRAL data we added in the model a normalization constant C. A value of C different from 
one would indicate a flux variation between the soft and hard-X ray energy domain. The derived 
values of C are well clustered around 1, with an average value of 1.45 and a standard deviation of 
0.76. For two sources only, namely IGR J16024-6107 and IGR J23308+7120, C is significantly 
larger than 1 (3.3^ and 3.1^', respectively), suggesting flux variability between the soft and 
hard X-ray observations. To support this conclusion we note that in the INTEGRAL/IBIS catalogue 
(Bird et al. 2010) both have a bursticity index flagged with yej^, and both show flux variation in 
the 20-40 keV band larger than a factor of three, in their ligh-curves. In addition, in the case of 
IGR J16024-6107 the high value of R would suggest that this object is peculiar in several aspects. 
By excluding these sources from the sample, the average value of C (standard deviation) becomes 
lower, 1.19 (0.31). As a result, we are confident that the value of R and E c we found are realistic 
and not due to variability. 

For the purpose of this work we do not examine the source variability in detail as we are mainly 
interested in explaining the average properties of the sample. However, to retain the completeness 
of our sample, we chose to include also the two variable sources (IGR J 16024-6 107 and IGR 
J23308+7120) in our analysis. 



The bursticity of a source is defined as the ratio of the maximum significance on any timescale, compared to the significance defined for the 
whole dataset. A flag of yes (Y) indicates a bursticity greater than 1. See Bird et al. 2010 for more details 
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Figure 4. Le/f panel: Distribution of the reflection fraction for all sources in the sample. Upper limits are shown whit obliques lines while detections 
are horizontal lines. Right panel: Relative reflection as a function of the photon index, filled boxes represent IGR sources, open circles are known 
AGN already known. 
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Figure 5. Left panel: Fe Ka line EW as a function of the absorbing column density N/f . The dashed line represents the EW values as obtained 
from transmission through a uniform slab and the solid line represents the function EW(N//)=EWoexp(<x/r e Nn), see text for details. Right panel: 
Reflection fraction R as a function of the scattering fraction f scan , filled boxes represent IGR sources, open circles are AGN already known. 



5.4 The Fe K line 

With the sole exception of IGR J17513-2011, a neutral Fe line is detected in all the sources of 
the sample (marginally significant for IGR J00040+7020, 95%). The value of the equivalent width 
(EW), calculated with respect to the total direct +reflected continuum^, ranges between 50-2000 eV 
(or 500 eV excluding the Compton thick sources). The average EW value considering all sources 
is 430 eV (with a somewhat high standard deviation of 564 eV), while, excluding the Compton 
thick objects, the average value becomes 196 eV (cr=176 eV). 



4 We calculated the EW as the ratio between the integrated photon flux of the Fe line and total (direct+reflected continuum flux at the energy of 
the line). 
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In the left panel of Fig.[5]we show the Fe Ka EW as a function of the absorbing column density. 
Below N//=10 23 5 cm -2 the EW trend is almost flat, with a best fit constant value of EW =120+16 
eV; this value is in a very good agreement with that found in other X-ray sample of Compton thin 
type 2 (Fukazawa et al. 2011) as well as type 1 AGN (Bianchi et al. 2009) observed with Suzaku, 
and XMM-Newton, respectively. However, above N#=10 235 cnr 2 the behaviour is different, and 
a clear correlation between N H and EW is evident, as expected from reflection from a Compton 
thick torus (Awaki et al. 1991). 

The question now is: what is the origin of the observed Fe K line in our sample? The dashed 
line in Fig. [5] represents the EW values as obtained from transmission through a uniform slab of 
medium with solar abundances subtending An to a continuum source of photon index Y = 2.0 
(Leahy & Creighton 1993). Large EW above 10 23 5 ~ 24 cm" 2 could be originated in a Compton 
thick torus, also responsible for the observed absorption. Since a typical torus with N#=10 24 cm -2 
should produce Fe(EW) =600-700 eV (Ghisellini et al. 1994) for solar abundances, to reproduce 
our values of about 1 keV, we have to consider a Compton thick medium with abundances higher 
than solar. On the other hand, the large EW observed in Compton thin AGN (N H below 10 23 5 cm -2 ) 
are far above the dashed line, suggesting that the line is due to reflection from a material different 
from the absorber. The evidence that the same Fe Ka EW (~120 eV) is found for Compton thin 
and type 1 AGN, suggests that the origin of this component is in the same medium, e.g. the BLR. 
The solid line in Fig. \5\ represents the function EW(N#)=EW exp(cr f >./V/ / ), with cr fe being the 
photoelectric cross-section at 6.4 keV and EW =120 eV. This trend well reproduces the Fe EW 
below 10 23 5 cm~ 2 , and represents a scenario in which the Fe line is generated by reflection off 
the BLR seen along an unobscured line-of-sight, while the underlying continuum is obscured by 
matter with a column density N H ', the total effect will be a larger value of the equivalent width 
increasing N#. In this regard, it is worth noting that in 4 sources of our sample (IGR J07565- 
4139, IGR J10404-4625, IGR J12026-5349, MCG-5-23-16) a broad Fe line is required to fit the 
broadband data. The intrinsic width of these broad lines is abroad ~0.3 keV corresponding to a 
velocity of v=15000 km s" 1 , that is completely compatible with an origin within the BLR. This 
scenario is supported by the recent observations of rapid variation of the absorbing column density 
in both Compton thin and thick AGN. In fact, a fraction of Compton thick (thin) objects do not 
intercept the torus along the line of sight, but are caught when a Compton thick (thin) cloud located 
in the BLR is passing in front of the source. Such sources may change their status in subsequent 
observations, once the thin or thick cloud has passed, explaining some of the so-called changing- 
look objects (the most clear case being NGC 1365, Risaliti et al. 2007). 
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Figure 6. The X-ray Baldwin effect, Iwasawa-Tanigichi effect. The EW of the Fe K line as a function of the 2-10 keV luminosity corrected for 
absorption (let panel) and the 20-10 keV luminosity (right panel). 



However, we stress that a Compton thin gas, also covering An solid angle to the source, cannot 
create a Fe EW of 50-1 50 eV (Yaqoob etal. 2001) as observed in our sample below N#=10 235 cnr 2 ; 
then, a Compton thick medium like a torus (not intercepting the line of sight in Compton thin or 
type 1 AGN), should exist. 

Our conclusion is that in heavily obscured Seyferts (N H above 10 23 5_24 cm" 2 ), the Fe line is 
produced in a torus also responsible for the absorption; for Compton thin sources (N H below 
10 23 5 cm" 2 ), as well as for type 1 AGN, a high fraction of the Fe K line (but not the total line pho- 
tons) has to be produced in a gas located closer to the black hole than to the torus. This region can 
be associated with the BLR that could also be responsible for the observed moderate absorption. 

We also investigated the X-ray Baldwin/'Twasawa-Taniguchi" effect (Baldwin 1977, Iwasawa 
& Taniguchi, 1993) indicating the existence of an anticorrelation between the 2-10 keV luminos- 
ity corrected for absorption and the EW of the Fe K line. This effect has been claimed in a sample 
of radio-quiet type 1 AGN (Bianchi et al. 2007, Page et al. 2004) in the form of EW aL™- 03 , 
with Lx,44 being the 2-10 keV luminosity in unit of 10 44 erg s _1 . This effect has not been revealed 
in a sample of type 2 AGN observed BeppoSAX data (Dadina 2008). However, the physical ex- 
planation for this phenomenon is still unclear and several possibilities are generally invoked in the 
literature, including a change of the ionizing continuum and gas metallicity with luminosity (Ko- 
rista et al. 1998) or a luminosity-dependent covering factor and ionization parameter of the BLR 
(Mushotzky & Ferland 1984). On the other hand, some studies claim that the primary physical 
parameter which drives the Baldwin effect could be the accretion rate, rather than luminosity. 

In the left panel of Fig. [6] we plot the EW of the Fe line as a function of 2-10 keV luminosity 
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corrected for absorption (L corr (2-10 keV)). A linear fit to the data of the form log(EW)=5.3+4.8 
- (0.07+0.1 l)log(L corr (2-10 keV)), indicates no evidence of a correlation (r=0.109). The X-ray 
Baldwin effect is not confirmed in our sample; although, we note that using the average EW in 
each luminosity bin, the correlation increases to r=0.2. In the right panel of Fig. [6l we plot the 
EW of the Fe line as a function of the hard X-ray luminosity in the 20-100 keV energy range. The 
correlation coefficient increases in this case to r=0.3, still indicating a lack of correlation. 

It is worth noting that the luminosity range we are exploring is rather narrow (log(L)=41.5- 
44.5), and this could make it difficult to measure any correlation as observed by Bianchi et al. 
(2009) over a wider luminosity range. 

5.5 Soft X-ray emission 

All sources in our sample show evidence for an unabsorbed component at energies less than 2 
keV. This soft component is likely produced by electron scattering of the intrinsic continuum 
and is observed in most Seyfert 2 galaxies (Antonucci 1993, Matt et al. 1997). To parametrize 
the fraction of the intrinsic flux scattered into the line of sight we used the f sc value. As shown in 
Table [2l our results suggest a value of f sc of a few per cent in all sources (with the exception of IGR 
J14515-5542 for which f sc =30% and IGR J16351-5806 with f sc =9%, see following discussion in 
this section). This value is in good agreement with that found in different samples of obscured 
AGN, observed with XMM-Newton and Chandra (Bianchi et al. 2006). 

Recently, Eguchi and collaborators (2009) discussed the evidence for an anticorrelation be- 
tween the amount of relative reflection R and the scattered fraction observed below 2 keV in a 
small sample of obscured AGN observed with Swift/BAJ and Suzaku. In particular the authors 
claimed the existence of a new population of type 2 AGN, characterized by high R (>0.8) and low 
scattered fraction f sc (< 0.5%) compared with the classical ones characterized by low R and high 
fs C . We investigated this possibility within our sample and in Fig.[5](right panel) we plot the value 
of relative reflection R with respect to the f sc . In our sample only one hard X-ray selected source, 
IGR J 10404-4625, is characterized by high R and low f sr ; nevertheless a linear correlation Pearson 
test produces a very poor correlation coefficient (r=0.09) indicating a lack of separation between 
two different populations. The anticorrelation might be interpreted as due to a degeneracy between 
different parameters in the fitting procedure (R and power-law normalization) as also pointed out 
by Comastri and collaborators (2010) who analysed a sample of obscured AGN observed with 
Suzaku. 
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The very soft X-ray spectra of six IGR sources (namely 07565-4139, 1009-4250, 14515-5542, 
16351-5806, 20286+2544, 23524+5842) clearly show more complexities than allowed by our 
BLM, showing either emission lines below 2 keV in addition to the simple scattered power-law 
component (see Figs. [9l CD} or a soft power-law slope significantly steeper than that of the 
hard X-ray continuum. Guainazzi & Bianchi (2007) found that a photoionized plasma (probably 
associated with the NLR) is an excellent description of the soft X-ray spectra of a large sample 
of obscured Seyfert 2 observed with the high resolution RGS gratings on board XMM-Newton. 
However, the available counting statistics and the limited CCD energy resolution of our present 
data are not sufficient to adequately constrain the physical properties of the emitting plasma. The 
short XMM-Newton exposures available, also make the RGS data for these four sources useless. 
We therefore adopted a simplified approach adding a thermal plasma component to the spectra of 
the four sources characterised by "soft-excess" emission. The temperature of this medium is in 
the range between 0.1-1 keV, while the total integrated soft X-ray luminosity of this component 
ranges from 10 40 erg s _1 to 1.4 x 10 41 erg s _1 . This excess emission could be associated to the 
photoionised NLR, unresolved sources and/or diffuse emission in the host galaxy (Guainazzi et al. 
2005). 



6 TYPE 1 VS TYPE 2 ACCRETION PROPERTIES 

Sample selection is the most important issue in testing the predictions of the UM scheme and 
several results inconsistent with the scheme can be explained as due to the biases inherent in the 
samples (Antonucci 2002). Several improved and well defined samples are attempting to test the 
validity and limitations of Seyfert unification (Singh et al. 2001, Ricci et al. 2011; Brightman & 
Nandra 201 1), however, issues related to the sample selection are still open. 

A broadband study of 41 type 1 AGN extracted from the hard X-ray selected complete sample 
(Malizia et al. 2009) has been presented by Molina et al. (2008). Following the UM, and having 
corrected the 2-10 keV luminosity for the absorption effect (L^To^y)' we would expect to find the 
same accretion properties (photon index, bolometric luminosity, Eddington ratio) between type 1 
and type 2 AGN belonging to a homogeneous sample. To investigate this issue we compared the 
results obtained for our absorbed sample and the unobscured one by Molina et al. (2008). 

From the samples of obscured and unobscured AGN, we collected all available published mass 
measurements obtained with reverberation mapping techniques, maser, X-ray variability or M-cr 
relation. A sub-sample of 33 (17 type 1 + 16 type 2) AGN with known masses has been obtained, 
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Table 3. Sub-sample of type 1+ type 2 AGN with measured Masses. 



src 


Mass 


reference 




(Mo) 




type 2 


NGC 788 


3.2e8 


Whang & Zhang 2007 


NGC1068 


10e7 


Whang & Zhang 2007 


NGC1142 


2.3 e9 


Winter et al. 2009 


MKN3 


3.0e8 


Whang & Zhang 2007 


MCG-5-23-16 


5e7 


Wandel & Mushotzky 1986 


NGC3281 


4.2e8 


Winter et al. 2009 


NGC4388 


8.5e6 


Kuo etal. 2011 


NGC4507 


1.8e8 


Nicastro et al. 2003 


NGC4945 


1.0e6 


Greenhill et al. 1997 


CenA 


2.0e8 


Marconi et al. 2004 


NGC5252 


9.5e8 


Capetti et al. 2005 


circinus 


1.7e6 


Marconi et al. 2004 


NGC6300 


2.e6 


Awaki et al. 2005 


ESO103-G35 


7.17e6 


Hayashida et al. 1998 


CygA 


2.1e9 


Marconi et al. 2004 


NGC7172 


8e7 


Awaki et al 2006 


Type 1 


mcg8-ll-ll 


1.5e7 


Bian & Zhao 2003 


MKN6 


1.7e8 


Winter et al. 2009 


NGC3783 


l.le7 


Peterson et al. 2004 


NGC4151 


4.57e7 


Peterson et al. 2004 


NGC4593 


9.8e6 


Peterson et al. 2004 


ESO511-G03 


4.6e8 


Winter et al. 2009 


MCG6-30-15 


4,5e6 


Bian & Zhao 2003 


NGC6814 


1.4e8 


Winter et al. 2009 


IC4329A 


2.2e8 


Peterson et al. 2004 


4C74.26 


le9 


Winter et al. 2009 


3C390 


3.9e8 


Peterson et al. 2004 


MR2251-178 


5.8e8 


Brunner et al. 1997 


IGR J07597-3842 


2e8 


Masetti et al. 2006 


IGR J12415-5750 


9.2e7 


Masetti et al. 2009 


IGRJ16558-5203 


7.8e7 


Masetti et al. 2006 


IGRJ16482-3036 


1.4e8 


Masetti et al. 2006b 


2E1853.7+1534 


1.4e8 


Masetti et al. 2006c 
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44 45 
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Figure 7. Left panel. Luminosity distribution (corrected for absorption) in (2-10 keV) for the sub-sample of type 1 and type 2 AGN with measured 
masses. Type 1 are shown as horizontal large spaced line, type 2 are diagonal small spaced lines. Right panel, the same distribution for the 20-100 
keV luminosity L 2 o-ioofev- 
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and reported in Table [3] In the case of unobscured AGN, Molina et al. (2008) showed that this 
sample can be considered as representative of the total one. For the obscured AGN in this work, 
the sub-sample we obtained is mostly representative of the known population, in fact we do not 
have any mass estimates for the new type 2 INTEGRAL AGN, and the sub-sample is composed 
by all the known sources in our original sample (excluding LEDA170194 and IC4518A). How- 
ever, we have previously shown (see Sect. 15.11 15.31 and 15.51) that IGR and the known objects are 
characterized by the same intrinsic properties. Therefore we are confident that the sub-sample can 
be representative of the whole one. We used this hard X-ray selected sample to draw a statisti- 
cal comparison between type 1 and type 2 Seyfert of (1) the absorption corrected 2-10 keV and 
20-100 keV luminosity, L corr (2-10 keV) and L(20-100 keV); (2) the bolometric luminosity, L^; 
and (3) the Eddington ratio, L foo //L ££W , with the main goal of testing the prediction of the UM 
for AGN. In Fig. [7] we show a comparison between the unabsorbed luminosities in 2-10 keV (left 
panel) and 20-100 keV (right panel) energy band. Horizontal large spaced lines and diagonal small 
spaced lines represent type 1 and type 2 class, respectively. The distribution of the 2-10 keV lu- 
minosity has averaged value (standard deviations) of 43.52 (0.71) and 43.06 (0.72) for type 1 and 
type 2, respectively. Using the Kolmogorov-Smirnov test we do not find a significant difference 
in the distributions, having a probability P=0. Ill to be drawn from the same parent population. 
This result is also confirmed in the higher, unbiased for absorption, energy range (20-100 keV) 
with averaged values (standard deviations) that are 43.82 (0.71) and 43.17 (0.70) for type 1 and 
type 2 respectively, with a KS test probability of 0.137. To be conservative and fully convinced to 
have properly corrected the 2-10 keV luminosities for absorption, we excluded the Compton thick 
sources from this analysis and, even in this case, the KS test probability confirms the similarity of 
the two populations (P=0.165). 

For all sources in the sub-sample we calculated the bolometric luminosity derived from the 2- 
10 keV luminosities, adopting the luminosity-dependent bolometric correction given in Marconi 
et al. (2004). In Fig. [H left panel,we show the distribution of hboi for type 1 and type 2 objects 
(same symbols used in Fig. [7]). The h bo i average value (standrad deviation) is 44.91 (0.91) and 
44.33 (0.90) for type 1 and type 2, respectively. A KS test gives a probability P=0.111 that the 
distributions belong to the same parent population, suggesting, even in this case, the similarity 
between the two classes of AGN. In the right panel of Fig. [8] we plot the photon index of the 
primary power-law versus the Eddington ratio for our sub-sample. The average value of the ratio 
~L bo i/L edd (with its standard deviation) is 0.17 (0.20) and 0.32 (0.72) for type 1 and type 2 AGN, 
respectively. Although the dispersion is huge, we note that, as for the luminosity, no evidence of 
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Figure 8. ie/i panel: Bolometric luminosity distribution for the sub-sample of type 1 and type 2 AGN with measured masses. Type 1 are shown as 
horizontal large spaced line, type 2 are diagonal small spaced lines. Right panel: Photon index as a function of the ratio l^tolf^edd for the sub-sample 
of type 1 and type 2 AGN with measured masses. Open triangles are type 1 while filled circles represent type 2 AGN. 



separation in the accretion rate distribution is evident between absorbed and unabsorbed AGN 
(KP probability P=0.292). These results concerning luminosities is also reflected in the spectral 
properties. In fact a comparison between the primary photon index in type 1 and type 2 class does 
not show any evidence for a separation (see right panel in Figure [8]) being the average value of Y 
equal to 1.80 (cr=0.18) and 1.75 (<x=0.26) for type 1 and type 2 respectively. A KS test in this case 
gives a probability P=0.862. 

These same conclusions have recently been reached by Singh and collaborators (201 1) through 
the X-ray analysis of an optical selected sample of type 1 and type 2 Seyfert. The value of the 
~L bo i/L edd ratio in our sub-sample of type 1 and type 2 AGN range in 0.001 - 3, extending to higher 
values with respect to a optical selected sample (Singh et al. 201 1). 

We are aware that this result could be affected by two issues. The first one concerns the mass 
estimation and its error, that can be as large as 10-20 per cent. It is still debated how accurate the 
[OIII] FWHM is as a proxy of the stellar dispersion <x and, especially for absorbed AGN, there 
are several inconsistencies between the different techniques used to measure the BH mass. For 
the second caveat we assumed equal bolometric correction for both type 1 and type 2 objects, 
however a different bolometric correction for absorbed and unabsorbed AGN (Barger et al. 2005) 
might affect our conclusions on the bolometric luminosity distributions. 

Finally, we also stress that, using the whole hard X-ray selected sample of 41 type 1+33 
type 2 from Molina et al. (2008) and the present work, we find that L corr (2-10 keV) and L(20- 
100 keV) have a small probability (P=0.05 and 0.02, respectively) to be drawn from the same 
parent population. This effect is likely to be due to a bias in the sample selection; in fact, due to 
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the less effective area that limits the sensitivity to only bright sources, samples of Seyfert AGN 
based on INTEGRAL, and Swift/BAT surveys have shown to contain a relatively large number of 
high luminosity and less absorbed Seyferts (Tueller et al. 2008; Treister et al. 2009; Beckmann 
et al. 2009). In fact, taking into account our total type 1 + type 2 (41+33) sample, we find that 
the average value of the redshift is 0.05 and 0.03 for type 1 and type 2 respectively, with type 1 
extending to higher redshift (higher luminosity) than type 2. 

Based on the model (Nicastro 2000) in which the BLRs are formed by accretion disc insta- 
bilities occurring in proximity to the critical radius at which the disk changes from gas pressure 
dominated to radiation pressure dominated, Nicastro and collaborators (2003) suggested that the 
absence or presence of hidden BLRs (HBLRs) is regulated by the ratio between the X-ray lumi- 
nosity and the Eddington luminosity, which, in the accretion power scenario, is a measure of the 
rate at which matter falls onto the central supermassive black hole. Following this model, a natural 
transition between type 1 and type 2 AGN should then occur, depending on the accretion rate. Ac- 
cording to Nicastro (2003), a value of the Eddington ratio of 0.001, i.e. the lowest value measured 
in our sub-sample of type 1+type 2 AGN (see Fig. [H right panel), represents the threshold value 
above which we expect to find HBLR (for obscured objects) or BLR (for the unobscured ones). 
Optical observations confirm this trend in type 1 AGN, while spectropolarimetric observations in 
a fraction of our sub-sample of type 2 AGN has to be performed to confirm the proposed scenario. 

7 CONCLUSION 

In this paper we have discussed the broadband X-ray properties of a sample of 33 absorbed hard 
X-ray selected AGNs. We selected a baseline model to reproduce the 0.3-100 keV data consisting 
of an absorbed primary emission with a high energy cutoff and its scattered fraction below 2-3 
keV, plus the Compton reflection features (Compton hump and Fe line emission). The broadband 
data coupled with the good energy resolution at the iron line energy range, allowed us to get a very 
good characterization of the primary continuum and of all the spectral features. 
The main results of this study can be summarized as follows: 

• We measure a high-energy rollover below 150 keV in 30% of the sample (10/33) and estimate 
a lower limit for the other sources is well below 300 keV. Our analysis suggests that the high energy 
cut-off is ubiquitous in the X-ray spectra of type 2 AGN. This evidence puts strong constraints on 
XRB synthesis model. 

• The hard X-ray selection allowed us to investigate the N H distribution in a sample much 
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less severely biased by absorption with respect to those selected in 2-10 keV band. When com- 
pared with 2-10 keV selected samples, our objects are characterized by larger column density; 
this evidence confirms the capability of hard X-ray observations for selecting obscured objects. In 
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addition, the high S/N ratio of our wide band data allowed us to isolate the Compton thick sources 
in a diagnostic plot (softness ratio F o „(2-10 keV)/F(20-100 keV) vs N H ). 

• The value of the Fe Ka EW as a function of the absorbing column density and the broad 
profile of the line in 4 sources, support the fact that (at least) a fraction of the Fe line has to be 



24 A. De Rosa et al. 




1 10 100 

Energy (keV) 

Figure 11. Continued 



produced in a Compton thin gas, likely associated to the BLR. This gas intercepts the line of 
sight and should also be responsible for the moderate absorption observed in the mildly absorbed 
sources (~10 23 cm -2 ). It is worth noting that we still need a Compton thick torus (not intercepting 
the line of sight) in moderately absorbed AGN to (1) produce the reflection continuum observed 
with INTEGRAL above 10 keV and (2) to contribute to the production of the Fe line emission. To 
confirm this scenario, broadband time-resolved variability study are needed, to directly observe 
the N H and Fe EW variations. These observational campaigns have been performed, so far, in just 
a few sources like NGC1365, NGC7582, i.e. the so called "changing look" AGN (Risaliti et al. 
2009 and reference therein, Bianchi et al. 2009). 

• A comparison with a hard X-ray selected sample of type 1 AGN shows that obscured and 
unobscured Seyfert are characterized by the same nuclear and accretion properties, i.e. luminosity 
(this result holds in both 2-10 keV and 20-100 keV bands), bolometic luminosity, Eddington ratio 
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and primary photon index. This evidence suggests that type 1 and type 2 objects should share the 
same central core. 

• The Eddington ratio would suggest, according to the model of Nicastro (2000), that all our 
type 2 AGN hide a BLR. Spectropolarimetric measurements are necessary to definitely establish 
the presence of HBLR clouds in the INTEGRAL AGN complete sample. 

More accurate broadband data up to 80 keV will became available in the near future with 
NuStar (Harris et al. 2010) to better investigate these issues. 
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